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A closed-form analytical technique has been developed to screen orbital average 
heating variations as a function of beta angle, altitude, surface area, and surface 
optical properties. Using planetary view factor equations for surfaces parallel-to 
and normal-to the local vertical, a cylindrical umbral shadow approximation, and a 
simplified albedo flux model, heating rate equations are formulated and then 
integrated to obtain orbital average heating. The results are compared to detailed 
analytical predictions using Monte Carlo integration and an assessment of error is 
presented. 

Nomenclature 

A = flat plate surface area 

a = albedo factor 

a = solar absorptivity 

(J = beta angle 

e = infrared emissivity 

h = circular orbit altitude above planet surface 

Q = heating rate 

<1 planet = planetary infrared flux 

q solar = solar flux 

r e = planet radius 

r = spacecraft radius vector 

s = solar vector 

0 - orbit angle measured from orbit noon 

VF = view factor to planet surface 

Sub- and Superscripts 

AVG = orbital average 

Z = zenith (space) facing surface 

N = nadir (planet) facing surface 

F = forward (velocity vector) facing surface 

A = aft facing surface 

P = port facing surface 

S = starboard facing surface 

|| = surface normal parallel to the local vertical 

J. = surface normal perpendicular to the local vertical 
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Introduction 

Characterization of the on-orbit thermal environment involves orbital parameters such as 
altitude and beta angle, environmental constants (solar, albedo, and planetary infrared 
constants), and spacecraft surface optical properties ( a and e ). Environmental heating 
calculations assessing the variation in one or more of these parameters has been 
relegated to the use of main frame computer analysis codes such as the Thermal 
Radiation Analyzer System (TRASYS) or Thermal Synthesizer System (TSS). Rickman 
and Ortiz (Reference 1) developed the Thermal Interactive Mission Evaluation System, 
which allowed for rapid calculation of simplified on-orbit thermal environments and 
parametric analyses. Such programs served as screening tools and their use 
demonstrated a significant reduction in the number of more detailed analyses that would 
have been required had the screening tool not been available. 

A useful screening calculation is presented in this work and can be expressed as a 
closed form solution given some constraining assumptions. Nevertheless, the described 
method allows for calculation of detailed environmental heating profiles as a function of 
orbit angle in addition to integrated orbital average heating rates for individual flux 
components. 

Assumptions 

This method described here is applicable given the following assumptions: 

a. the central body has a constant infrared flux emission over its entire surface; 

b. the albedo factor, a , is assumed constant over the entire surface; 

c. the spacecraft is orbiting at a low altitude (i.e., h«r e ) and in a circular orbit; 

d. the overall spacecraft time constant is on the order of the orbit period; 

e. the spacecraft is oriented in a fixed local-vertical local-horizontal attitude with 

surfaces facing in the principal directions as discussed below. 

Assumption (a) is required to simplify the integration of the planetary heating flux and is 
applicable for Earth-orbiting spacecraft. 

For assumption (b), while the abledo factor is constant, the albedo flux varies as the 
cosine of the orbit angle. 

Assumption (c) constrains the analysis to low altitudes in order to allow the simplifying 
assumption of a cylindrical umbral shadow and a locally constant albedo heating flux. 
These concepts will be developed in greater detail below. 

Assumption (d) validates the use of orbit average heating as a screening tool, however, 
the expressions established for heating as a function of orbit angle are useful for 
screening instantaneous heating to surfaces. 
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Assumption (e) greatly simplifies the view factor analysis since only two types of 
planetary view factors need be calculated. This assumption, however, limits the 
applicability of this technique to surfaces facing in the principal directions (i.e., forward-, 
aft-, starboard-, port-, zenith- and nadir-facing surfaces). 

General Expressions 

Consider a spacecraft in a low altitude, orbit about a planet as presented in Fig. 1. Due 
to precession of the orbit plane and the planet’s motion about the sun, the angle, p, 
between the solar vector, s , and its projection onto the orbit plane will vary as a function 
of time. 


zp 



Fig. 1 Problem geometry and p definition. 

As p varies, so too will the fraction of the orbit in which the spacecraft spends in the 
planet’s umbral shadow. At Earth’s distance from the sun, the umbral shadow cone 
stretches over 800,000 miles into space. Hence, the shadow cone is nearly cylindrical in 
shape close to Earth and deviates by only about 0.26° from the assumed cylindrical 
shadow. With this simplification, an expression relating the terminator entry and exit 
angles, Gentry and 0 Exn -, respectively, with p is easily obtained. This simplifying 
assumption does not provide for the transition through the penumbra. At low altitudes, 
however, this does not pose a problem as the time spent in the penumbra is measured in 
seconds compared to an orbit period on the order of 1.5 hours. 

As depicted in Fig. 2, the position of the spacecraft, r , can be expressed as a function of 
the altitude above the planet, h, the radius of the planet, r e , the angle from orbit noon.G, 
and p: 


r = ( r e + h)cos0 cos J3i +(r e +h)s\n9j + (r e +/i)cos0sin J3k 

The projection of this vector onto the YpZp-plane is: 

r' = (r e + h)smdj' + ( r e + h) cos# sin fik 


(1) 


( 2 ) 
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And the magnitude of the projection is simply: 


|r '| = ( r e + h)^j sin 2 # + cos 2 # sin 2 


( 3 ) 


Vp 



Fig. 2 Geometry of the simplified eclipse problem. 

Examination of the geometry reveals that the onset of the umbral shadowing occurs 
when the magnitude of ?' is less than r e , or: 


sin# < 


r. 


' COS 2 P 


K r e +h 


-sin 2 p 


(4) 


Since sin 0 is defined between -n/2 and +n/2, we note that 0 ENT ry = n-0 and 0 Exn - = 7t+0. 

The heating incident on an orbiting flat plate varies as a function of the orbit angle, 0, and 
is given by: 


A PLATE 
X TOTAL 


(»hQ 
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(5) 


For six plates, arranged as shown in Fig. 2, the total absorbed heating as a function of 0 
is given by: 


QtOTAl{@)~ £&*«(*)+ y . Q ALBEDO iP ) + ^jQ PLANET iP) (®) 

m=Z,N,F,A.P,S m=Z,N,F,A,P,S m=Z,N,F,A,P,S 

where the superscripts Z, N, F, A, P, and S refer to the zenith-, nadir-, forward-, aft-, 
port-, and starboard-facing plates, respectively. 
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View factor calculations are required in order to determine the amount of planetary 
infrared and albedo fluxes impinging on the analytical surfaces. Reference 2 provides a 
thorough discussion of planetary- and albedo-radiation view factor calculations for a 
planet-oriented and arbitrarily-oriented flat panel. However, the treatment presented 
herein is restricted to surfaces oriented normal-to or parallel-to the nadir direction with 
the further restriction that h « r e . Reference 3 provides the general derivation of the 
planetary view factor for a flat plate surface and has been adapted, here, for the special 
cases of plates whose surface normals are parallel (nadir-facing) to and perpendicular 
(forward-, aft-, port-, and starboard-facing) to the local vertical vector, respectively: 


(W0„ - 


\ r '+ h , 




n- 2sin 


-i 


\ r e+ h S 


-sin 


2sin" 


1- 


r. ^ 


■Y\l 


\ r e +h J 
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Note, also, that since the zenith facing surface cannot view the planet: 


(VF) 


ZENITH 


= 0 


(7) 

( 8 ) 

(9) 


Solar Heating 


The orbital average absorbed solar heating is obtained by direct integration of the 
instantaneous solar heating equations over the applicable angles and dividing by the 
angular span of an entire orbit (27i). For orbits that pass into the planet’s umbral shadow, 
i.e., when: 






\ 2 


\ r e+ h J 


-sin 2 p 


>0 


the integrals are: 


{(2 solar )avg 


= ^QME^z cos l 

2 n J -' r/2 


( 10 ) 


( 11 ) 


{Q solar )avg 
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(13) 
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For p > 0: 


Iq y _ <isoLAR a P A P s i n ft f ENTRY dQ 

SOLAR JAVG ^ L 

Z7T EXIT 

(15) 

O 

II 

£ 

1 

§ 

(16) 

For p < 0: 


{Q SOLAR )aVG “0 

(17) 

(q Y - ~^SOLAR a S^S si* 1 /? f ENTRY dQ 

SOLAR JAVG 2 ^- J?/ v// 

(18) 

and for orbits that experience no eclipse, the limits of integration 0 ENTRY 

and 0 EX | T given 


above may be replaced by +71 and - 71 , respectively. 

A Simplified Albedo Model 

A true albedo model would account for variations in the albedo constant as a function of 
solar incidence angle as well as variations in the view factor to a surface with varying flux 
intensity. The effect of the latter of these two factors is most evident near the ground 
terminators where the slope of the cosine reduction in solar flux on the planet’s surface is 
greatest. For low altitude orbits, however, it is assumed that only a small portion of the 
planet is visible and that any local variation in albedo flux intensity is small. As an 
example, at 220 nm altitude, the distance to the horizon is approximately 1200 nm. At a 
P of n/2 (i.e, 90°), the spacecraft is orbiting along the ground terminator with a view 
subtending 0.35 radians in all directions. Assuming a pure cosine reduction in the 
albedo flux, the maximum intensity of the albedo flux, at the extent of the view, is 
approximately 34% of that at the subsolar point. Since half of the plate is viewing a 
darkened planet, an upper limit of about half that value, or 17%, of the subsolar albedo 
flux is assumed in the worst case. Since the visible portion of the planet’s surface is not 
uniformly lit (with the albedo flux going to zero at the ground terminator) a more realistic 
estimate of the albedo error in this extreme case is closer to 10%. It should be noted 
that analyses performed at lower values of p will exhibit less error since little time is 
spent at or near the ground terminators. At the subsolar point, the variation over the 
fraction of the visible illuminated surface is, at most, about 6%. 

With the inherent error in hand, then, the simplified orbital average albedo heating can 
be expressed as the sum of the integrals over the illuminated portion of the orbit (-7t/2_< 0 
< +ti/ 2) divided by angular span of the entire orbit (27i): 

({2 ALBEDO istl'G = ® 

(Qalbsoo L = fcos m (20) 

2 n • L?r/2 
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(h Y - a ( VF )^ SOLAR® F A F C0S l V'lfun 

albedo Ja vg ~ 2 n 1^/2 
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)-rtH 
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Planetary Infrared Heating 


The planetary infrared heating is easily calculated due to the assumption of uniform flux 
emission over the entire planet. Due to the constant planetary infrared heating 
assumption, the orbital average planetary infrared heating is not a function of0 and is 
given by: 

(QpLANEt)aVG - 0 

(25) 

(QpLANEt)aVG planet € n^n 

(26) 

{Q PLANET )a vg - (^)l 9 PLANET £ F 

(27) 

{Q PLANET L VG = (^O-L Q PLANET £ A ^ A 

(28) 

{G planet) avg = (^)i? planet 8 pAp 

(29) 

{Q PLANET ^AVG ~ (^)-L # PLANET £ S 

(30) 

Results 


The equations presented above are integrated and the results are 

presented on a 


component-by-component basis. 

For the solar heating where an eclipse occurs: 
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{Q solar) AVG ~ 0 ^ SOLAR 01 F^F COS P 
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For p > 0: 


For (3 < 0: 
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(35) 

(36) 

(37) 

(38) 


and for the condition when no eclipse occurs, the integral for the zenith surface remains 
the same and the other surfaces become: 


(39) 
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For p > 0: 


For p < 0: 


{Q solar )avg ~ QsoLAR a pAp sin P 
\Q SOLAR) AVG 

{Q SOLAR )aVG = 0 

{Q SOLAR JaVG = SOLAR a s ^S Sm P 


(42) 

(43) 

(44) 

(45) 


(33) 
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For the albedo heating: 
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And, finally, for the planetary infrared heating: 
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In order to test the method, a sample case was developed and compared to detailed 
Monte Carlo predicted heating rates calculated using the Thermal Synthesizer System 
Heatrate application (Reference 4). The simple geometry of a unit box (1 ft x 1 ft x 1 ft) 
with optical properties for all box sides set to a=1 and e=1 was selected so as not to 
improperly bias the results away from a given viewing direction and heating component 
with smaller inherent error. The box was assumed to orbit the Earth at an altitude of 
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220 nm inclined 67° with respect to the equator and would experience the entire range of 
possible p angles. The solar flux, albedo factor, and planetary infrared flux were 
assumed to be 443.7 Btu/hr ft 2 , 0.3, and 77.0 Btu/hr ft 2 , respectively; all well within the 
accepted ranges (Reference 5). Detailed Monte Carlo analyses were performed for a 
range of p using 10000 rays per surface for each heating component. A comparison of 
the orbital average heating resulting from both methods is presented in Fig. 3 with error 
presented in Table 1. 



Beta Angle (Degrees) 

Fig. 3 Comparison of simplified screening technique and detailed Monte Carlo 
predicted average heating for an orbiting box at 220 nm altitude. 


Beta (deg) 

Closed-Form 
Solution (Btu/hr) 

Detailed Solution 
(Btu/hr) 

% Error [ 100 x (Detailed - 
Closed)/Detailed ] 

-90 

600.1 

607.4 

1.2 

-80 

706.4 

707.8 

0.2 

-71 

787.9 

788.6 

0.1 

-70 

784.2 

786.5 

0.3 

-60 

691.9 

692.3 

0.1 

-40 

682.8 

682.6 

0.0 

-20 

653.9 

654.2 

0.1 

0 

581.7 

581.9 

0.0 

20 

653.9 

654.2 

0.1 

40 

682.8 

682.6 

0.0 

60 

691.9 

692.3 

0.1 

70 

784.2 

786.5 

0.3 

71 

787.9 

788.6 

0.1 

80 

706.4 

707.8 

0.2 

90 

600.1 

607.4 

1.2 

Table 1 1 

fabular results and 

a comparison of the detailed and 


closed-form solution error. 
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Examination of the results indicates that the simplified method produces exceptionally 
good results when compared with the detailed predictions. Errors on the order of 1% are 
observed at p = 90° with much lower error at lower values of p. These errors are well 
within the statistical error observed in the detailed Monte Carlo analysis results. 

Potential Applications 

The method appears to be well suited for screening where approximate solutions will 
suffice. Establishing the worst-case hot and cold orbits prior to a more detailed analysis 
has the potential to significantly reduce the magnitude of the overall thermal analysis 
cycle. If implemented in spreadsheet form, engineers can rapidly investigate the effects 
of optical property degradation, altitude and beta angle variations. Expressions for 
heating as a function of orbit angle are useful for studies of instantaneous orbit heating. 
Such an algorithm may also have application as a screening tool in a larger analysis 
system and would promote the development of an “intelligent system” to pre-select 
analysis cases. 

Concluding Remarks 

A simplified, closed-form method for screening orbital heating variations as a function of 
the p, altitude, surface area, and optical properties has been developed and compared to 
detailed analytical predictions. The method is easily implemented in computerized 
spreadsheet programs and provides a rapid assessment of orbital heating trends, 
identifies maximum and minimum overall heating conditions, and allows exploration of 
the sensitivity of the spacecraft absorbed heating characteristics through variation of 
design parameters such as optical properties and areas. Implementation of this method 
may prove useful as a screening tool for spacecraft thermal engineers prior to the 
initiation of detailed orbital heating analyses. 
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